Analog of the spin-orbit induced anomalous Hall effect with quantized radiation 
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We demonstrate how the term describing interaction between a single two-level atom and two 
cavity field modes may attain a Rashba form. As an outcome, cavity QED provides a testbed for 
studies of phenomena reminiscent of the spin-orbit induced anomalous Hall effect. The effective 
magnetic field, deriving from the non-Abelian gauge potentials rendered by the Rashba coupling, 
induces a transverse force acting on the phase space distributions. Thereby, the phase space dis- 
tributions build up a transverse motion manifesting itself in spiral trajectories, rather than circular 
ones obtained for a zero magnetic field as one would acquire for the corresponding Abelian gauge 
potentials. Utilizing realistic experimental parameters, the phenomenon is numerically verified, 
ascertain that it should be realizable with current techniques. 
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Progress in manufacturing and a deeper understanding 
of electronic devices have caused an increased interest in 
the field of spintronics [l|. For example, the spin-orbit 
coupling occurring in various materials such as semicon- 
ductors provides a handle for the electron spin. In fact, 
the spin-orbit coupling is a main ingredient for both the 
spin Hall effect [2[ and its companion, the anomalous Hall 
effect (AHE) [|. This makes these effects attractive from 
a practical point of view, and not only as phenomena on 
their own right. However, despite its long history, the 
underlying physics of especially the latter, the AHE, has 
remained a controversial subject [|| . 

Already E. Hall discovered that the Hall effect is 
greatly enhanced in ferromagnetic compared to non- 
magnetic materials This became known as the AHE. 
It is only recently, with knowledge of geometrical phases, 
that a more complete picture of the AHE has been ob- 
tained d. The processes contributing to the AHE can be 
divided into three categories: (i) intrinsic (deriving from 
the band structure, and consequently from the spin-orbit 
coupling), (ii) skew scattering (originating from distorted 
scattering of impurities), and (iii) side-jump (often iden- 
tified as the contribution making up for the difference be- 
tween the total AHE conductivity and the intrinsic plus 
skew conductivities). While the first has a topological 
origin, the last two contributions are due to impurities in 
the material. Numerous measurements of the AHE con- 
ductivity have been presented |6j, but identifying, and 
even controlling, the different mechanisms underlying the 
AHE current remains an extremely non-trivial task. In 
the realm of spintronics, it would therefore be desirable 
to come up with alternative systems free from impuri- 
ties such that the AHE current derives solely from the 
spin-orbit coupling, which could, for example, lead to 
accomplishing effective Datta-Das transistors [7j. Over- 
coming the obstacle of impurities with solid state devices 



seems very difficult, and one should preferably look in 
other directions. 

Two possible candidates, both known to exhibit Hall 
characteristics, are either rotating Bose-Einstein conden- 
sates [H or cold atoms in optical lattices @. Due to 
the purity and high controlled versatility of these sys- 
tems, great interest has been paid to them during the 
recent past The AHE has been predicted to oc- 

cur on the p-band of cold fermions in an optical lattice 
An optional system considers cold atoms coupled to 

in which the spin Hall 
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spatially varying light fields 

effect has been predicted fl~3l ] . Regardless of the devel- 
opment of novel experimental techniques in cooling and 
trapping cold atom gases, presently non of the different 
Hall effects nor non-Abelian properties have been realized 
utilizing cold atoms. 

Yet another system endowed from impurities and pos- 
sessing long coherence time scales is cavity quantum 
electrodynamics (QED); a single atom interacting with 
an isolated set of cavity modes [ijj. During the last 
two decades, cavity QED has successfully demonstrated 
among others, entanglement generation [l5j . the quan- 
tum measurement problem and the quantum-classical 
transition [16j . and verification of the graininess of the 
quantized electromagnetic field [l7[- Moreover, recent ex- 
periments have realized coherent coupling of single quan- 
tum dots 18 1 or Bose-Einstein condensates [l|| to a cav- 
ity mode, paving the way for the possibility of reaching 
a super-strong coupling regime of cavity QED. 

Notwithstanding the many advantages provided by 
cavity QED models, relatively few schemes employing 
them as quantum simulators have been suggested. In 
Ref. [2(| it was presented how effective non-Abelian 
gauge potentials arise in bimodal cavity QED models. 
Similar cavity QED systems had also been shown to 
mimic Jahn- Teller models frequently appearing in molec- 
ular and condensed matter physics [2l| . A proposal how 
to realize magnetic monopoles by the means of cavity 
QED systems was outlined in [22|]. A model similar in 
structure, the cold trapped ion system, was analyzed 
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in [231 ] and it was put forward how it may simulate rela- 
tivistic particles (see also [13] )■ in this paper we demon- 
strate the appearance of transverse phase space currents, 
evocative of intrinsic AHEs, in a bimodal cavity QED 
system. 

We consider a high-Q cavity containing a two-level 
atom (quantum-dot) dipole-interacting with two degen- 
erate but orthogonal cavity modes. The system Hamil- 
tonian reads 



H = H f + H a + B 1 



where 



H f = huj (Va + 6 + S) , H a =J2 E i\3){j\, 
Hi = d- P(x). 
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Here, a) and w (a and b) represent the creation (anni- 
hilation) operators for the two modes, lu is their com- 
mon mode frequency, Ej the internal atomic energy of 
state \j), d is the dipole transition moment, and P(x) is 
the electric field. In terms of the internal atomic states, 
the components of the dipole moment operator becomes 
d a = — e| l)ct{2| — e|2)a(l|, where e is the electron charge 
and a = x,y, z. In the dipole approximation we set 
x = 0, and the field is written as 



E — E\£\i (a — d^ 1 ) + E^E-ii [b—P 



(3) 



where Ej is the polarization vector for mode j and Ej 
its corresponding field amplitude. A deeper insight of 
the system characteristics is obtained by expressing the 
fields in their quadrature operators 



P = -= (a - a f 
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obeying the canonical commutation relations; [Xk,Pk'\ — 
iSkk' ■ For brevity we introduce the Pauli operators a x — 
|1)(2| + |2)(1|, * y = -i|l)<2|+i|2)<l|, a, = - |2)<2|, 
and set the zero energy such that E\ = Ml/ 2 and E2 — 
— Ml/2. Furthermore, we choose g\ = d-E\E\\[2 = —ga x , 
and assume gi = d ■ E2E2V2 = —go y with g real. The 
resulting Hamiltonian is the E x e Jahn- Teller one [2l| 
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first introduced in the molecular physics community [25| . 
Via a simple 7r-rotation around the a x axis, the atom- 
field interaction term attains a Rashba form [2(| . In the 
absence of A 12 -terms, the same type of Hamiltonian is 
frequently used to describe a two dimensional gas of free 



spin-orbit coupled electrons. Adding the parabolic po- 
tential renders, on the other hand, a stereotype model 
representing spin-orbit coupled quantum dots [27| . Such 
confinement naturally restricts the particle motion, but 
nonetheless, as will be shown spin-orbit induced trans- 
verse currents still play a crucial role for the system dy- 
namics. 

The adiabatic potential surfaces (APS) of 
He x*, defined as Vj( fj.ft) = Mj (Pf + P 2 2 ) /2 ± 
hy/tt 2 /4: + g 2 {Pf +P$), are envisaged in Fig. [U These 
have a polar s ym metry and possess a conical intersection 
at the origin [25|]. The cr z -term split the degeneracy at 
the origin, and whenever g < s/cjfl the sombrero shape 
of the lower APS is lost and instead a global minimum 
at the origin is attained. Conical intersections appearing 
in momentum space are frequently referred to as Dirac 
points due to their linear dispersions. 




Figure 1: The APSs of the Exe Hamiltonian ([5]). The conical 
intersection is located at the origin, Pi = P2 =0. 

The Hamiltonian is conveniently written on the form 

„ I (r, -- -I,-"' 

Hexc 




where 



(Ai,A 2 ) = — (a x ,a y ) 



U) 



(6) 



(7) 



The quantities Ak and $ are effective vector and scalar 
potentials respectively. The gauge invariance follows di- 
rectly from considering their response to unitary trans- 
formations; 

d 

A -> UUX, t)AU(X, t) - P f (X, t) — U(X, t), 

(8) 

$ W(x, t)$u(x, t) - ihu f (x, t)—u(x, t). 

Pointed out in Ref. (28|, the intrinsic AHE is a direct 
result of a non-zero Berry phase originating from a vector 
potential (also called Mead-Berry curvature). One way 
to illustrate the appearance of the Hall effect is by the 
generalized spin-dependent Lorenz force analogous to the 
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work by Wong on color-charged classical particles moving 
in non-Abelian fields [2{|. As an outcome of the gauge 
potentials, there is an associating magnetic field 

Bi = ^s ijk F kl , F kl = d k Ai - diA k - i[A k ,A{\. (9) 

The first two terms of F^i are identically zero in our 
model. On the other hand, the last term, deriving from 
the non-Abelian structure of the system [29], is indeed 
non-zero. The effective magnetic field, being propor- 
tional to <5>, induces a state (spin) dependent force; the 
magnetic field is either in the positive or negative z- 
direction depending on the internal atomic state |1) or 
1 2). Due to the opposite direction of the force on the in- 
ternal states, the system may act as a spin-filter (Datta- 
Das transistor) which explains its interest for spintronics. 
Here, however, the intrinsic transverse AHE current does 
not take place in real space as in solid state devices, but 
in a phase space representation of the field states. 

For typical atom cavity QED parameters it is legiti- 
mate to apply the rotating wave approximation, which 
would cast the Hamiltonian ([5]) into a bimodal Jaynes- 
Cummings one [30[. On the other hand, for circuit cavity 
QED employing solid state quantum dots/SQUIDs in- 
stead of true atoms, imposing the rotating wave approx- 
imation is usually not justified. In fact, in the parameter 
regimes where this approximation is applicable one finds 
that the transverse current rendered by the spin-orbit 
coupling is rather weak, and spontaneous emission of the 
atom as well as cavity decay will thereby become signif- 
icant within the time needed to build up a noticeable 
transverse current. 
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Figure 2: (Color online) The upper plot (a) displays the aver- 
age radius of the phase space distributions as function of time. 
Black curves shows the ideal situation without losses, while for 
the light blue curve losses have been taken into account. The 
two lower plots present the same time evolution but of (-Xi) 
vs. (X2) without (b) and with (c) losses. The system param- 
eters are P w = 0, P20 = 0, X10 = 10, X20 = 0, Q/2ty = 6.9 
GHz, uj/2-k = 5.7 GHz, g/2n = 105 MHz, 7/2tt = 1.9 MHz, 
and k/2-k = 250 kHz. 

Normally, the transverse AHE current is driven by an 



electric field having the effect of tilting the 2D Air- 
plane. In this cavity QED setting, a comparable tilt 
is generated by external driving of the cavity. Pump- 
ing with an amplitude 77 of say mode 1, implies a term 

Hp Ump = hi] (o-i + o^ij = h\[2r\X\ in the Hamilto- 
nian ©. This induces a position shift of oscillator 1. 
An alternative way of picturing the pumping is to dis- 
place the initial field state of mode 1 by the operator 
D(a) — e aa i~ a Sl , with a = r}\2/u. In our numerical 
simulations we will consider this alternative viewpoint 
by assuming the initial state of mode 1 to be displace 
by a, whereby we omit any pumping term Hp Ump in our 
numerical calculations. 

In a solid, impurities prevents the electrons to accel- 
erate to infinity causing an equilibrium situation. It is 
known that a constantly accelerating situation results in 
no net transverse current (3lj . The present model lacks 
impurities, but on the other hand the confining potential 
of Eq. © restricts the motion in the AiA2-plane and the 
spin-orbit induced AHE is thereby still apparent in the 
phase space motion. 

To demonstrate the parallel of an intrinsic AHE in the 
present system we solve the Schrodinger equation em- 
ploying the split-operator wave packet method. In par- 
ticular we utilize the experimental parameters of the cir- 
cuit QED experiment presented in Ref. [l8j], and more- 
over assume the atom (quantum-dot) to be initially in its 
ground state |1), mode 2 in vacuum (0)2, and mode 1 in 
a coherent state |aio)i (displaced vacuum); 

= V ip i° Xl e-< Xl - Xl °) 2 / 2 (10) 

Here, X w and P10 are respectively the initial position 
and momentum of the coherent state, related to the am- 
plitude as aio = (A10 — iPw)/y2. By choosing P w = 
and Xio > 0, the combined field wave packet starts out 
with a zero momentum at (Xi,X2) = (A"io,0). Apart 
from the potential force causing the wave packet to slide 
down the harmonic potential, in addition it experiences 
the effective transverse Lorentz force originating from 
the magnetic field (O. For an Abelian gauge potential 
(<7i oc 172), the Lorentz force vanishes and the field wave 
packet would remain localized along the Ai-axis through- 
out the propagation. In the present model, however, the 
oscillatory motion will depart from the Ai-axis inducing 
a non-zero field in mode 2. As time progresses, popula- 
tion is transferred from mode 1 to mode 2 until finally 
mode 1 is in vacuum. Consequently, in phase space, mode 
1 will follow an inward spiral trajectory while mode 2 fol- 
lows an outward trajectory until the population is com- 
pletely swapped between the two modes after some time 
t s . This is analogues to the intrinsic AHE; considering 
polar coordinates r and ip instead, one has that the radial 
motion along r induces a transverse rotational current 
in the (^-direction. For longer time periods, the process 
starts over again with the population of the two modes 
interchanged. The spiral motion is depicted by the black 
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curves in Fig. [2] (a) showing the averages (Xi) 2 + (Pi) 2 
for initial values X 2 q = Pw = P20 = 0, X w = 10, and 
the atom initially in |1). Figure [2] (b) displays the corre- 
sponding time evolution of (A"i) versus (A2). At about 
80 ns (for the utilized parameters) , the oscillations are al- 
most entirely along the AVaxis indicating the swapping 
of population between the two modes. 

So far we considered a closed system neglecting sys- 
tem losses. From Fig. [2] it is clear that the two modes 
perform a very large number of individual oscillations 
within the time t s . Thus, decoherence of both the cav- 
ity fields and the atom may become important over such 
time periods. To investigate the effect of losses, we con- 
sider the effective non-hermition Hamiltonian i? e // = 

Hexc ~ i K (j^a + — i7|2)(2|. Here k and 7 are the 
photon decay rate and the atomic spontaneous emission 
rate respectively. The results corresponding to the black 
curves of Fig. [2j but where experimental values of k and 
7 have been considered, are shown as light blue curves in 
Fig. [5J The AHE clearly survives despite the relatively 
long interaction time, only a slight drop in field inten- 
sity, reminiscent of momentum relaxation, is found. It 
is understood that this analysis of the effect of an envi- 
ronment does not correctly capture decay of coherence of 
the system state. This, however, should not be a problem 
in the present situation as we are here only interested in 
the field amplitude of the two modes and not coherent 
superpositions nor the precise values of the quadrature 
operators. 



Summarizing, in this work we have shown how trans- 
verse phase space currents, reminiscent of AHEs, natu- 
rally occur in bimodal cavity QED systems. The phe- 
nomenon has been explained in terms of an effective 
non-Abelian gauge field and its corresponding state- 
dependent Lorentz force acting upon the combined phase 
space distributions. In particular, the effective gauge 
potentials induces phase space trajectories having spiral 
shapes; starting in coherent states, one field will follow an 
inward spiral trajectory and the other an outward spiral 
trajectory causing in general a swapping of population 
between the two modes. Utilizing experimental parame- 
ters, we showed that the intrinsic AHE should be visible 
in current experiments even in the presence of losses. Fi- 
nally, we point out that bimodal cavity experiments have 
been performed [331 ] . and moreover that, apart from field 
intensities, even the field quadrature operators ((4]) are 
easily measured. Indeed, the ENS group of S. Haroche re- 
cently presented experimental results were the full phase 
space distribution of a cavity mode was assessed pjij . 
It should also be mentioned that the general idea pre- 
sented in this work is not restricted to a two-level atom. 
It equally well applies to other atom-field configurations, 
such as for example a three-level A-atom coupled to two 
cavity modes. 
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sions, and acknowledge support from the MEC program 
(FIS2005-04627). 
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